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Abstract 

In a search for signatures of physics processes beyond the Standard Model, various 
eeqq vector contact-interaction hypotheses have been tested using the high-Q^, 
deep inelastic neutral-current e+p scattering data collected with the ZEUS detector 
at HERA. The data correspond to an integrated luminosity of 47.7 pb~^ of e~^p 
interactions at 300 GeV center-of-mass energy. No significant evidence of a contact- 
interaction signal has been found. Limits at the 95% confidence level are set on the 
contact-interaction amplitudes. The effective mass scales A corresponding to these 
limits range from 1.7 TeV to 5TeV for the contact-interaction scenarios considered. 
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1 Introduction 



The HERA ep collider has extended the kinematic range available for the study of 
deep inelastic scattering (DIS) by two orders of magnitude to values of Q"^ up to about 
50000 GeV^, where is the negative square of the four-momentum transfer between the 
lepton and proton. Measurements in this domain allow new searches for physics processes 
beyond the Standard Model (SM) at characteristic mass scales in the TeV range. A 
wide class of such hypothesized new interactions would modify the differential DIS cross 
sections in a way which can be parameterized by effective four-fermion contact interac- 
tions (CI) which couple electrons to quarks. This analysis was stimulated in part by an 
excess of events over the SM expectation for > 20000 GeV^ reported recently by the 
ZEUS and HI collaborations, for which CI scenarios have been suggested as possible 
explanations (see e.g. 



In a recent publication [0], the Born cross sections at > 400 GeV^ extracted from 
47.7 pb^^ of ZEUS neutral-current (NC) e'^p DIS data collected during the years 1994- 
1997 have been compared to the SM predictions primarily derived from measurements at 
lower Q^, extrapolated to the HERA kinematic regime. The agreement is generally good. 
The only discrepancy is due to the high-Q^ excess of events in the data taken before 1997, 
which has not been corroborated by the 1997 data but is still present in the combined 
sample with reduced significance. The present paper presents a comparison of the full 
data sample to SM predictions modified by various hypothetical eeqq CI scenarios. Limits 
on the CI strength and on the effective mass scale A are determined for the different CI 
types. 

Limits on eeqq CI parameters have been reported previously by the HI collaboration [§, 
by the LEP collaborations (ALEPH 0, DELPHI 0, L3 in|,|ll, OPAL ll|), and by 



the Tevatron experiments CDF [|T^] and D0 |17] 



This paper is organized as follows: after a synopsis of the relevant theoretical aspects in 
Sect. H, the experimental setup, the event selection and reconstruction procedures, and 
the Monte Carlo (MC) simulation are discussed in Sect. |^. The CI analysis methods 
are presented in Sect. and the results are summarized in Sect. A discussion of the 
statistical issues related to the limit setting and the information needed to combine the 
results of this analysis with those of other experiments are given in the Appendix. 



2 Contact— Interaction Scenarios 



A broad range of hypothesized non-SM processes at mass scales beyond the HERA center- 
of-mass energy, ^/s = 300 GeV, can be approximated in their low-energy limit by eeqq 



1 



contact interactions (Fig. ||), analogous to the effective four-fermion interaction describing 
the weak force at low energies fl^. Examples include the exchange of heavy objects with 
mass M 2> a/s such as leptoquarks or vector bosons and the exchange of common 
constituents between the lepton and quark in compositeness models ||21],^ . Note that the 
CI approach is an effective theory which is not renormalizable and is only asymptotically 
valid in the low-energy limit. 

In the presence of eeqq CIs which couple to a specific quark flavor (g), the SM Lagrangian 
£sM receives the following additional terms p0|-p3|: 



'nt{^LeR){qLqR) + Vs'i^L^R)i^RlL) + h.c. scalar 
+ vlLi^LYeL){qL'yf^qL) + vlR{eLl^eL){qRl^qR) 
+ TlRL{^Rl^eR){qL-if,qL) + VRR{eRl^eR){qR-ff,qR) vector 



(1) 



+ r]Ti^L(y^''eR){qL(y^^qR) + h.c. 



tensor 



where the subscripts L and R denote the left- and right-handed helicity projections of 
the fermion fields, g is the overall coupling, and A is the effective mass scale. Since g 
and A always enter in the combination g"^ / h?, we adopt the convention g^ = An so that 
CI strengths are determined by the effective mass scale, A. The overall sign of the CI 
Lagrangian is denoted by e in eq. (|ID. Note that e = +1 and e = — 1 represent separate CI 
scenarios which are related to different underlying physics processes. The t] coefficients 
determine the relative size and sign of the individual terms. Only the vector terms are 
considered in this study since strong limits beyond the HERA sensitivity have already 
been placed on the scalar and tensor terms (see [^^Q and references therein). In the 
following, only CI scenarios will be discussed for which each of the t]"^^ (m, n = L,R) is 
either zero or ±1. 

The relevant kinematic variables for this analysis are Q^, x and y, which are defined in 
the usual way in terms of the four-momenta of the incoming positron {k), the incoming 
proton (P), and the scattered positron {k') as = —q^ = —{k — k')'^, x = Q'^/{2q ■ P), 
and y = {q ■ P)/{k ■ P). The leading-order, neutral-current ep SM cross section is given 

by 

^^^(x, Q^) - [(1 + (1 - vr) ^ (1 - (1 - vr) ^^r] , (2) 

Fr{x,Q') = Yl A(Q')- N^,Q') + a:g(x,Q^)] , (3) 

q=d,u,s,c,b 

xF^ix^Q^) = Yl B^{Q')-[xq{x,Q')-xq{x,Q')] , (4) 

q=d,u,s,c,b 

where q{x, Q"^) and q{x, Q^) are the parton distribution functions for quarks and anti- 
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quarks, a is the fine structure constant, and Ag and Bg are defined as follows: 



2 (5) 

B,{Q')= [{Vg-){A-g) - {Vg-){A- 



Terms resulting from CIs can be included in eqs. (|2|-^ by replacing the SM coefficient 



functions V^'^ and A^'^ with 



= Qg- {V, ± ae) Vg xz + ^(^^L + vLr) 



eQ2 



f 



1 

4sin2^^cos2^^Q2 + jvf| ■ 

In eq. (^), the subscript m is L or R; the plus (minus) sign in the definitions of and 
A^ is for m = L {m = R). The coefficients and a/ are the SM vector and axial- 
vector coupling constants of an electron (/ = e) or quark (/ = ?); Qf and Tj denote the 
fermion's charge and third component of the weak isospin; Mz and 6^ are the Z mass 
and the Weinberg angle. In the limit A oo, the coefficient functions 
eq. (PD reduce to their SM forms. 



As can be seen from eqs. (H-y), the effect of a CI on the NC DIS cross section depends 
on the specific scenario. In general, two kinds of additional terms are produced. One 
kind is proportional to and enhances the cross section at high Q^. The second is 
proportional to 1/A^ and is caused by interference with the SM amplitude, which can 
either enhance or suppress the cross section at intermediate Q^. The predicted ratio 
(SM+CI)/SM of (Pa /{dQ"^ dx) depends on at fixed a;, but also on x at fixed due to 
the different y dependences of the coefficient functions multiplying the F2 and -F3 terms 
in eq. @j. Note that CIs induce modifications of the SM cross section for all x and Q^, 
in contrast e.g. to the direct production of an eq resonance in the s-channel. 

For ep scattering at HERA, the contribution of second- and third-generation quarks to CI 
cross-section modifications is suppressed by the respective parton distribution functions 
in the proton. For the present analysis, flavor symmetry, 

Vmn Vmn Vmn and Vmn Vmn C"^) 

is assumed unless explicitly stated otherwise. The CI limits reported here are only weakly 



3 



sensitive to this assumption.^ Contributions from the top quark content of the proton are 
almost completely suppressed due to the large top mass and are neglected in this analysis. 

Using the relations in eq. (0), there are eight independent vector terms in eq. (|l]), which 
lead to a large list of possible CI scenarios. To reduce this list, we consider the following: 



• Recent measurements of parity-violating transition amplitudes in cesium atoms |^ 
imply very restrictive constraints on CIs |^,^,^. These limits are avoided by parity- 
conserving CI scenarios, i.e. if 

vIl + vIr - v'rl - = . (8) 

Conforming to this constraint, in particular, excludes CIs of purely chiral type, i.e. 
those for which the rj'^^ are non-zero only for one combination of m and n. 

• SU{2)l invariance requires 77^^ = 77^^, 0- Terms violating this relation are considered 
only for u quarks, which dominate the high-a; cross section at HERA, and hence show 
the largest CI-SM interference effects for a given A. A CI signal from this source 
could therefore manifest itself in ep collisions while avoiding strong SU (2)-breaking 
effects e.g. at LEP. 

Based on these considerations, the 30 specific CI scenarios listed in Table are explored in 
this paper. Note that each line in this table represents two scenarios, one for e = +1 and 
one for e = —1 (denoted as VV"*", VV~ etc.). All scenarios respect eq. (|), and all scenarios 
except U2, U4 and U6 obey SU{2) symmetry. The S'[/(2)-conserving CI scenarios with 
Vll 7^ Vll (Ul and U3) would also induce an euqq CI signal in charged-current (CC) DIS, 
e+p — i> VX. We have not used the CC data sample to constrain further these scenarios. 

Several examples of modifications of the SM cross sections by CIs are illustrated in Fig. ^. 
The cross-section modifications for the X1-X6 and the corresponding U1-U6 scenarios 
are similar, demonstrating that the d quarks have little impact on the CI analysis. 



3 Experimental Setup and Data Samples 

This analysis uses the data samples, Monte Carlo simulation, event selection, kinematic 
reconstruction, and assessment of systematic effects used in the NC DIS analysis described 
in 0. The data were collected during the years 1994-1997 in e~^p collisions with beam 
energies Ee = 27.5 GeV and Ep = 820 GeV. The relevant aspects of the experimental 

^ A similar statement is true for the Tevatron CI limits from lepton-pair production, which also depend 
on the parton distributions in the proton. In contrast, the CI analyses at LEP are sensitive to the cross 
section a{e'^e^ —^ hadrons) and the resulting limits depend strongly on flavor symmetry assumptions. 
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setup, event selection, and reconstruction are summarized briefly below. More details can 
be found in 0. 



The ZEUS detector is described in detail elsewhere The main components used in 
the present analysis are the central tracking detector (CTD) [Q, positioned in a 1.43 T 
solenoidal magnetic field, and the compensating uranium-scintillator sampling calorimeter 
(CAL) 131|], subdivided into forward (FCAL), barrel (BCAL) and rear (RCAL) sections. 



Under test beam conditions, the CAL energy resolution is a{E)/E = 18%/a/-E' [GeV] 
for electrons and 35%/^^ E [GeV] for hadrons. A three-level trigger is used to select 
events online. The trigger decision is based mainly on energies deposited in the calorime- 
ter, specifically on the electromagnetic energy, on the total transverse energy, and onQ 
E — pz = Yli-^ii^ ~ cos 6i) (the sum running over all calorimeter energy deposits). For 
fully contained events, the expected value oi E — pz is given by 2Ef, = 55 GeV. Timing 
cuts are used to reject beam-gas interactions and cosmic rays. 

The luminosity is measured to a precision of 1.6% from the rate of energetic bremsstrah- 
lung photons produced in the process ep ep^ . 

The offline event reconstruction applies an algorithm to identify the scattered positron 
using the topology of its calorimeter signal and the tracking information. The mea- 
sured energies are corrected for energy loss in inactive material between the interaction 
point and the calorimeter, for calorimeter inhomogeneities, and for effects caused by redi- 
rected hadronic energy from interactions in material between the primary vertex and the 
calorimeter or by backsplash from the calorimeter (albedo). The kinematic variables for 
NC DIS candidate events are calculated from the scattering angle of the positron and from 
an angle representing the direction of the scattered quark. The latter is determined from 
the transverse momentum and the E — pz of all energy deposits except those assigned to 
the scattered positron. 

The appropriately corrected experimental quantities are used to make the offline event 
selection. The major criteria are |^: (i) the event vertex must be reconstructed from the 
tracking information, with |Zvtx| < 50 cm; (ii) an isolated scattered positron with energy 
E'^ > 10 GeV has to be identified; (iii) 38 GeV < E - pz < Q5 GeV; (iv) < 0.95, where 
He is the value of y as reconstructed from the measured energy and angle of the scattered 
e"*". The requirements (iii) and (iv) reject background events from photoproduction. 

Monte Carlo simulations are used to model the expected distributions of the kinematic 
variables x, y and and to estimate the rate of photoproduction background events. 
NC DIS events including radiative effects are simulated using the HERACLES 4.5.2 p7 



^ ZEUS uses a right-handed Cartesian coordinate system centered at the nominal interaction point, 
with the Z axis pointing in the proton beam direction. The polar angle is defined with respect to this 
system in the usual way. 
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program with the DJANGO 6.24 |39] interface to the hadronization programs. In HERA- 
CLES, corrections for initial- and final-state radiation, vertex and propagator corrections, 
and two-boson exchange are included. The underlying cross sections are calculated in 
next-to-leading order QCD using the CTEQ4D0 set of parton distribution functions 
(PDFs). The NC DIS hadronic final state is simulated using the color-dipole model of 
ARIADNE 4.08 and, as a systematic check, the meps option of lepto 6.5 @ for 
the QCD cascade. Both programs use the Lund string model of jetset 7.4 [^] for the 
hadronization. MC samples of photoproduction background events are produced using the 
HERWIG 5.8 pO] generator. All MC signal and background events are passed through the 



detector simulation based on geant |51[], incorporating the effects of the trigger. They 
are subsequently processed with the same reconstruction and analysis programs used for 
the data. All MC events are weighted to represent the same integrated luminosity as the 
experimental data. 

Good agreement is found in [0 both between the distributions of kinematic variables in 
data and MC, and between the measured differential cross sections da/dQ"^^ da/dx and 
dojdy and the respective SM predictions, with the possible exception of the two events 
at g2 > 35000 GeVl 



4 Analysis Method 

The CI analysis compares the distributions of the measured kinematic variables with the 
corresponding distributions from a MC simulation of events of the type e"^p e'^X, with 
the weight 



^(SM+CI) 



w 



d?c 



dx dQ- 



■(SM) 



(9) 

true x,y,Q^ 



applied to each reconstructed MC event in order to simulate the CI scenarios. The weight 
w is calculated as the ratio of leading-ordei]^ cross sections, evaluated at the "true" 
values of x, y and as determined from the four-momentum of the exchanged boson 
and the beam momenta. In cases where a photon with energy is radiated off the 
incoming positron (initial-state radiation), the e beam energy is reduced by the energy 
of the radiated photon. The reweighting procedure using eq. (|^) accounts correctly for 
correlations between the effects of a CI signal and the pattern of acceptance losses and 
migrations. 



^ The final versions of the CTEQ5 Q and MRST Q PDF sets became available only after completion 
of this analysis. 

^ Note that CIs are a non~renormalizable effective theory for which higher orders are not well-defined. 
Radiative corrections due to real photon emission are expected to cancel to a large extent in eq. (0). 
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The simulated background events from photoproduction are added to the selected NC-DIS 
MC data sets. The photoproduction contamination is highest at high y and is estimated 
to be less than 0.5% overall and below 3% in any of the bins used for the cross-section 
measurements in |^. 

For each of the CI scenarios, two statistical methods are used. Each incorporates a log- 
likelihood function^ 



where the pi are appropriately normalized probabilities which are derived from a com- 
parison of measured and simulated event distributions {i runs over individual events in 
method 1 and over bins of a histogram in method 2, see below). Note that the two CI 
scenarios corresponding to two sets of r^^^ values differing only in the overall sign e are 
combined into one log-likelihood function. A description of the data samples used for 
evaluating L(e/A^) is given in Table |^ for both methods. 

• Unbinned fit to the (x, ?/)-distribution: 

The available experimental information entering the analysis can be split into two 
parts, the shape of the (x, distribution, {dP N / dx dy) / N tot , and the total number of 
events, Atot- The latter is related to the total cross section (Ttot(e/A^) in the kinematic 
region under study by iVtot = Ctot ■ ■ o'tot, where etot is the average acceptance and £ 
denotes the integrated luminosity. 

In the first method, an unbinned log-likelihood technique is applied to calculate 
Li(e/A^) from the individual kinematic event coordinates This method only 

makes use of the shape of the (x, y) distribution. The sum in eq. (0) runs over all 
events in the selected data sample. The MC events are appropriately reweighted to 
simulate a CI scenario with strength e/A^, as outlined in eq. (^). The probability den- 
sity p{x,y; e/ A'^) required to calculate pj(e/A^) = e/A^) is determined from 
the resulting density of MC events in x and y and is normalized to unity, thereby 
discarding the information on C ■ crtot(e/A^). The justification for this deliberate re- 
duction of experimental information is given a posteriori by the fact that crtot(e/A^) 
depends only weakly on A in the parameter space of interest: for A values larger than 
the 95% lower exclusion limits (see Sect. |), crtot(e/A^) deviates from the SM value, 
c"tot(0), by less than 2% for all scenarios except XI and X6, for which 2.2% and 2.6% 
are reached, respectively. This sensitivity is smaller than, or of the same order as, the 
1.6% systematic luminosity uncertainty quoted in Sect. |^ and is hence not significant. 

^ A discussion of a probabilistic interpretation of the log-likelihood function based on a Bayesian ap- 
proach can be found in the Appendix. 




(10) 



7 



Even though A^tot is not used in this method, it is important to note that data 
(-^tot = 13243 events observed) and SM prediction (13151 events expected) agree 
within the luminosity uncertainty and the statistical error. Furthermore, Figs. ^ and|^ 
demonstrate that the and x distributions agree in shape with the SM expectation. 
Therefore, a deterioration of the agreement of data and expectation with increasing 
|e/A^|, indicated by an increase in the log-likelihood function with respect to a min- 
imum close to e/A^ = 0, can be interpreted in terms of CI exclusion limits on e/A^ 
or A. 

The sensitivity of the results to systematic effects is studied by repeating the limit 
setting procedure (see below) for analysis parameters and selection requirements which 
are varied within admissible ranges. These systematic checks include those which were 
performed in the underlying cross-section analysis 0: 

— use of MC samples generated with the meps instead of the ARIADNE option (as 
described in Sect. H); 

— variations of trigger or reconstruction efficiencies and of experimental resolutions 
within their uncertainties by suitably reweighting the MC events; 

— modifications of the cuts and parameters used for event reconstruction and selec- 
tion; 

— variation of the calorimeter energy scales in the analysis of the data but not in 
that of the MC events. 

In addition, systematic uncertainties related to the CI fitting procedure are investi- 
gated by 

— use of MC samples which were generated with the following alternatives to the 
CTEQ4D PDF set: (i) with the PDF set MRS A [||] and (ii) using the resuhs of 
a recent NLO QCD fit [||] to the 1994 ZEUS H and HI H structure function 
data and to fixed-target data; 

— use of MC samples generated with the PDF sets CTEQ4A2 and CTEQ4A4 |4 



corresponding to values of as(M|) = 0.113 and 0.119 (instead of as(M|) = 0.116 
used for CTEQ4D); 

— changing the amount of photoproduction background in the MC sample by ±100%; 

— modifying details of the method used to infer p(x, y; e/A^) from the MC event 
distributions; 

— calculating the CI cross sections in the following alternative ways: 

— i> with different parton distribution functions, 

— > using NLO instead of LO QCD calculations and parton distributions, 

with the couplings restricted to first-generation quarks {j]mn—"^m.n—'^ran — O); 
— > with the couplings restricted to first- and second-generation quarks (?7m„=0). 

The procedure to determine limits on the CI parameters including the information 
from the systematic checks is discussed at the end of this section. 
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• Binned fit to the distribution: 

In the second method, L2{e/ A^) is determined from the distributions, using Poisson 
statistics for the numbers of events in each interval. Here, the sum in eq. ([To|) runs 
over all bins. For the calculation of L2, both the shape and the normalization of 
dN/dQ"^ are used. 

The systematic uncertainties are included in L2 using the assumptions that they are 
fully correlated between bins and that the probability densities for all uncertainties 
have Gaussian shapes. The effects taken into account are equivalent to those described 
above for the unbinned method and include in addition a 1.6% uncertainty on the 
integrated luminosity. 



Both methods have been shown to provide unbiased estimates of the CI strength when 
applied to MC samples. A few examples of comparisons of the resulting log-likelihood 
functions and L2 are shown in Fig. |^. Both functions agree with each other for most of 
the CI scenarios under study. However, for a few scenarios such as X2~ and X3~, rises 
faster than with decreasing A. This can be understood as a consequence of the fact that 
Li uses the full two-dimensional information of the (a;, y) distribution and is hence more 
sensitive to those scenarios which imply a marked x-dependence of the modification to 
the cross-section at fixed . The two-minimum structure of the log-likelihood functions 
seen in the AA and XI scenarios in Fig. ^ is characteristic for several CI scenarios for 
which the destructive SMxCI interference term cancels approximately the pure CIxCI 
term in a range of typically 0.1 TeV"^ < je/A^] < 0.2TeV~l MC studies indicate that the 
exact shape of 2 in the vicinity of these double-minima is dominated by the random 
pattern of statistical fiuctuations of the event distributions, but that the two-dimensional 
method has a higher probability than the one-dimensional method to assign a larger value 
of L to the "non-SM minimum" than to the "SM-minimum" (the AA case in Fig. ^ is 
typical). This is again understood as a consequence of the additional input information 
for the two-dimensional method. The normalization information, Nioi/{C ■ (Ttot(e/A^)), 
cannot distinguish between the minima since the difference of crtot(e/A^) between them is 
much less than the uncertainty of C 

The best estimates, Aq, for the different CI scenarios are given by the positions of the 
respective minima of Li 2(e/A^) with e = —1 and e = +1. The values of e/A^ resulting 
from the unbinned method are indicated in Fig. ^ Note that 1/A„ = is taken for all 
cases where Li_2(e/A^) rises monotonically for a given scenario, i.e. for e = +1 or e = — 1 
(see Fig. i. 

An analysis of the log-likelihood functions (referred to as the "L-analysis" in the fol- 
lowing) is usually employed to calculate confidence level intervals (i.e. limits) of e/A^. A 
discussion of some aspects related to the L-analysis approach can be found in the Ap- 
pendix, where polynomial parameterizations for Li(e/A^) are also provided, which may 
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prove useful for combining our results with those from other experiments. One problem 
with the L-analysis is that its results depend on the choice of the "canonical variable" 
in terms of which it is performed; for example the A limits are different if L is evaluated 
as a function of e/A instead of e/ h?. In order to avoid this ambiguity, MC experiments 
(MCE) are used, i.e. statistically independent MC data samples corresponding to the data 
luminosity. The presence of CIs in the MCEs is simulated by reweighting the events in 
the MCEs according to eq. @. For each MCE, the log-likelihood analysis is performed 
as a function of the assumed "true" value of A, Atrue, for each of the CI scenarios under 
study. The lower limit of A at 95% C.L. for a given CI scenario is determined as the value 
of Atrue at which 95% of the MCEs produce most likely values of |e/A^| larger than that 
found in the data. In the cases where l/A^ 7^ 0, the MCEs are also used to estimate the 
probability, psM, that a statistical fluctuation in an experiment with the SM cross section 
would produce a value of Aq smaller than that obtained from the data. Note that a high 
value of does not in itself signify that the SM prediction describes the data well, but 
indicates that the inclusion of the CI scenario under study does not significantly improve 
the agreement between data and prediction. 

For the two-dimensional fitting method, the above procedure is repeated for each sys- 
tematic check, using statistically independent MCE sets which reflect the corresponding 
modifications of the analysis. Each such MCE set consists of 500 MCEs. The resulting 
CI limits are scattered around the limits of the central analysis, with deviations in both 
directions being about equally frequent .0 The A limits deviate from their central values 
by typically less than 15%, though by as much as 30% in a few cases. The modification of 
the underlying SM cross section induced by a variation of as(M|) and by using different 
PDF sets (see above) causes variations of the A limits of typically a few percent and 25% 
maximally. 

Systematic effects are finally taken into account in the CI limit analysis by combining 
the MCE sets of all systematic checks and determining the values of Atrue for which 95% 
of all MCEs in the combined set produce most likely values of |e/A^| larger than that 
found in the data. This procedure is an approximation to averaging over the spectra 
of systematic effects, assuming that the different checks are uncorrelated and that the 
ranges of parameter variations (e.g. of the calorimeter energy scales or of q;s(M|)) reflect 
the actual uncertainties. The corresponding question defining a 95% C.L. limit is: "which 
value of Atrue causes deviations from the SM prediction which are larger than that observed 
in the experimental data in 95% of all ZEUS-type experiments exhibiting systematic 
differences according to the spectra determined in the analysis of systematic effects" .[] 



^ This implies that roughly 50% of all checks produce limits which are stronger than the central ones and 
is related to the fact that, to a very good approximation, A depends linearly on continuous parameters, 
like as, within their uncertainty intervals. 

^ This corresponds to assigning equal a priori probabilities to each of the tested variations and reflects 
the fact that by construction neither of them can be excluded or favored. 
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The A limits resulting from both log-likelihood methods agree to within 15% in all cases 
except for the scenarios AA"*", X2~ and X3~, for which the two-dimensional method has 
higher sensitivity and correspondingly yields significantly stronger limits. Therefore, the 
results of this method are presented in the following. 

5 Results 

The resulting SM probabilities (see Table ^ do not indicate significant amplitudes 
for any of the CI scenarios considered. Therefore, we report upper limits on 1/A^ and the 
corresponding lower limits on A. 

A selection of plots demonstrating the expected modifications of the Q^- and a;-distri- 
butions in the presence of CIs with strengths corresponding to the 95% CL. exclusion 
limits is shown in Figs. ^ and As mentioned in Sect. |1|, the data show an excess over 
the SM predictions at > 35000 GeV^. However, it is apparent that CIs cannot provide 
an improved description of this excess while simultaneously describing the data well at 
lower Q^, where data and SM expectation are in good agreement. These figures also 
confirm that the x-dependence of the (SM+CI)/SM cross-section ratio differs markedly 
between different CI scenarios and obviously contributes to the sensitivity of the CI fit, 
e.g. in the X3 case. This statement is generally true for all cases where the limits derived 
from the two analysis methods differ significantly. 

The lower limits on A (Ajj^^^) and the probabilities are summarized in Table § and are 
displayed in Fig. ^. In none of the cases does the SM probability fall below 16%. The A 
limits range from 1.7 TeV to 5 TeV. Those few cases with limits below 2.5 TeV correspond 
to log-likelihood functions having a broad minimum, either in the region with e = +1 
(XI, Ul) or with e = — 1 (X6, U6); these minima correspond to parameter combinations 
for which the pure CIxCI contribution and the CIxSM interference term approximately 
cancel in the HERA kinematic regime.^ 

Table ^ shows a comparison of the ZEUS CI results with corresponding limits reported 
recently by other experiments which study eeqq CIs in e+e" scattering at LEP (ALEPH 
0, L3 |12[, OPAL [0) or via Drell-Yan pair production in pp scattering (CDF |T6| , 
D0 [|l3)- The HI ||^ and DELPHI [|l^] collaborations report results only for purely chiral 
CIs which cannot be compared to the results of this paper. All limits shown in Table ^ 
have been derived assuming flavor symmetry (see eq. (0)), except the LEP limits for the 
U3 and U4 scenarios which are for first-generation quarks only. Limits for the X2, X5, 
Ul, U2, U5 and U6 scenarios are not included in Table ^because there exist no previously 

^ Note that the binning of the data in x and used for this presentation is irrelevant in the analysis 

of the limits provided in Table ^. 
^ Note that these cancellations happen at opposite e for e^p scattering. 
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published results. ZEUS and the other experiments are all sensitive to CIs at mass scales 
of a few TeV. The relative sensitivity to different CI scenarios depends on the CIxSM 
interference sign which is opposite in e~^p scattering on the one hand and in e~^e~ and 
pp scattering on the other.[^ Where available, the LEP limits often exceed the results of 
this paper, although it should be noted that this depends on the assumption of a flavor- 
symmetric CI structure. Limits for CIs which couple only to first-generation quarks would 
differ only by a small amount from those reported here for ep or pp scattering, but would 
be significantly weaker in the case of LEP. 

6 Conclusions 

We have searched for indications of eeqq contact interactions in 47.7 pb~^ of ZEUS high- 
e~^p neutral-current deep inelastic scattering data. The distributions of the kinematic 
variables in the data have been compared to predictions derived for 30 scenarios of vector 
contact interactions which differ in their helicity structure and quark-flavor dependence. 
In none of the cases has a significant indication of a contact interaction been found and 
95% C.L. upper limits on have been determined for each of these scenarios. The lower 
limits on A range between 1.7 TeV and 5 TeV and are found to be largely independent of 
the statistical method applied. 

The results exhibit a sensitivity to contact interactions similar to that recently reported 
by other experiments; in order to allow full use to be made of the available experimental 
data, the information needed to combine the results of this analysis with those from 
other sources is provided. Some of the limits reported here are the most restrictive yet 
published, and several of the contact-interaction scenarios have been studied in this paper 
for the first time. 
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The CI limits from CDF and the LEP experiments have been quoted here according to the sign 
convention used in the cited papers. 
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Appendix: The Log— Likelihood Functions 



In this Appendix we summarize some aspects of interpreting the log-hkehhood functions 
L(e/A^) using a Bayesian probabihstic approach. The results of the unbinned method 
described in Sect, ^jhave been employed here, but systematic effects have not been taken 
into account. For simplicity, we will denote the log-likelihood functions by L instead of 
Li in this Appendix. 

For ease of calculation, the functions L(e/A^) — Lmin have been parameterized as eighth- 
order polynomials in the region where L{e/ h?) — Lmin < 18, corresponding approximately 
to a ±6cr interval around the minimum of L. The polynomial coefficients are summarized 
in Table ^ The accuracy of the parameterizations is typically better than 0.1 units in 
L. Note that, neglecting systematic effects, these parameterizations allow one to combine 
the ZEUS results with those of other experiments by simply adding the L functions and 
repeating the analysis described below. 

The Bayesian approach starts from the relation 



where V symbolizes the experimental data, p{T>\t/K^) is the conditional probability to 
observe T) for a given value of e/A^, and Po{e/A'^) is the prior probability describing the 
knowledge about e/A^ before the experiment was conducted. The probability p(e/A^|D) 
assigned to e/A^ under the condition of having observed V is what we actually want to 
derive. 

In the following, we will identify 



with the normalization appropriately fixed to unity. In the simplest case of a Gaussian 
probability distribution, L(e/A^) is a parabola, and a^/Aa, the RMS width of p{V\e/A'^), 
corresponds to the width of the Gaussian. Even though some of the L(e/A^) functions of 
the CI analysis are not parabola-like, o'^/a2 is still well defined and can be interpreted as 
a measure of the experimental sensitivity to a given CI scenario. The values of o"e/A2 are 
summarized in Table |^. 

Usually, simple assumptions about pQ^e/Is?) are made in order to calculate p(e/A^|P), 
which only weakly depends on these assumptions provided that the width of po{e/h?') 
is much larger than a^/\2. We have calculated p(e/A^|X>) using a fiat prior probability 
restricted to either e > or e < 0. One-sided 95% C.L. limits in the Bayesian approach 



p{e/K^\V)^p{V\e/K^).p,{e/A^), 



(11) 



p(D|6/A2)ocexp(-L(e/A2)) , 



(12) 
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(^Hm b) have been determined by solving^ 

1/(A+ )2 

°- = 0.95 and 7' = 0.95; (13) 

/ rfeexp(-L(0) / rfeexp(-L(0) 



-oo 



For all CI scenarios, the results deviate by less than 20% from the limits resulting 
from the MCE method (see Sect. §). Note that systematic effects have not been considered 
for this cross check. 

An alternative way to present the results of this search analysis is obtained by considering 
the ratio of two equations of the type (0) for different values of e/A^, where one is taken 
as a reference value (chosen to be e/A^ = 0, i.e. corresponding to the SM). Rearranging 
the terms yields 

p{e/A'\V) / p{e/A' = 0\V) _ p{V\e/A') ^ 

Po(e/A2) / p^(,/A2 = 0) p(P|e/A2 = 0) ^ ^ ^ ' ^ ' 

where the double ratio on the left-hand side quantifies how the probability assigned 
to a given value of e/A^ changes due to the experimental data T>, with the reference 
point, e/A^ = 0, fixing the normalization. The function R{t/A^) has been discussed in 
detail elsewhere (see e.g. and references therein). The representation of eq. (|T^ is 
independent of the prior probability and can be used to combine the results of this analysis 
with those of other experiments by analyzing the product of corresponding R functions. 
Note that R does not involve integrations over e/A^ and is hence invariant with respect to 
the variable transformations mentioned in Sect. ^; in particular, R can be interpreted both 
as a function of e/A^ and as a function of A. By definition, R asymptotically approaches 
unity if e/A^ — or A oo, indicating the loss of experimental sensitivity as the CI 
strength vanishes. Regions where R is close to zero are excluded by P, whereas R > 1 
indicates "signal-type" regions where the experimental data are better described by a 
CI scenario than by the SM. Typical R values for significant deviations from the SM are 
expected to exceed unity by several orders of magnitude (cf. the discussion in [Q). Two 
representative examples of the functions -R(A) observed in the CI analysis are shown in 
Fig. 1^. In none of the scenarios from Table |1] does R exceed unity by more than 40%, 
corroborating the conclusion that there is no significant indication for the presence of 
CIs. The threshold-type region of R{A) indicates the position of the A limit. Indeed one 
obtains hmits similar to those reported in Sect, ^by solving the condition R{K^^ j^) = 0.05 



Similar approaches have been used by ALEPH |^ and CDF 
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V'rl 


1 nix 


Vll 


Vlr 


Vrl 


Vrr 


VV 




+1 


-Li 


+1 

1 -L 


+ 1 

1 -L 


-Li 

1 -L 


-Li 

1 -L 


-1-1 


AA 


+ 1 

1 -L 


— 1 


— 1 


+ 1 

1 -L 


+ 1 

1 -L 


-1 


— 1 


-Ll 

1 -L 


VA 


+ 1 


— 1 


-Li 


— 1 


+ 1 


— 1 


-Li 


— 1 


XI 
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—1 








-Li 


— 1 








X2 







-Li 
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1 -L 





-1-1 







-Li 


n 


n 


+1 


-Li 


n 


n 


-Li 

\ -L 


X4 


n 


-Li 


-Li 


n 


n 


-Li 


-Li 


n 


X5 





+ 1 





+1 





+ 1 





+1 


X6 








+1 


-1 








+1 


-1 


Ul 


+1 


-1 




















U2 


+1 





+1 

















U3 


+1 








+1 














U4 





+1 


+1 

















U5 





+1 





+1 














U6 








+1 


-1 















Table 1: The 30 scenarios for contact interactions considered in this paper. Each 
row of this table corresponds to two different CI scenarios for overall interference 
signs e = +1 and e = —1, respectively (see eq. in the text). 



quantity 


method 1 


method 2 


Qlin (GeV^) 


500 


400 


QL. (GeV^) 


90200 


51200 




0.04 




■^max 


0.95 




?/min 


0.04 




Umax 


0.95 




events 


13243 


37379 



Table 2: The kinematic regions used for the CI analysis in the two fitting methods. 
Note that the cuts on Q^, x, and y indicated in the table are applied in addition to 
the event selection criteria described in Sect. |^. 
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CI 


e = -1 


e = +1 


0"e/A2 


^lim 

(TeV) 


PSM 


(TeV) 




vv 


5.0 




4.7 


0.28 


0.021 


VA 


2.6 


0.25 


2.5 


0.25 


0.070 


AA 


3.7 


0.28 


2.6 




0.080 


XI 


2.8 


0.26 


1.8 




0.113 


X2 


3.1 




3.4 


0.28 


0.056 


X3 


2.8 




2.9 


0.37 


0.066 


X4 


4.3 




4.0 


0.26 


0.034 


X5 


3.3 




3.5 


0.28 


0.052 


X6 


1.7 


0.16 


2.8 


0.27 


0.105 


Ul 


2.6 


0.24 


2.0 




0.125 


U2 


3.9 




4.0 


0.38 


0.037 


U3 


3.5 




3.7 


0.48 


0.046 


U4 


4.8 




4.4 


0.32 


0.025 


U5 


4.2 




4.0 


0.36 


0.032 


U6 


1.8 




2.4 


0.24 


0.118 



Table 3: Values of the 95% C.L. limits (kf^J for all CI scenarios, as well as 
the SM probabilities, psm, for those scenarios with e/A^ ^ (see Sect. In the 
last column, the RMS width of the probability distribution pi oc exp(— L(e/A^)) is 
indicated (cf. Appendix). 
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CI 


A^^ (TeV) (95% C.L.) 


ZEUS'" 
this 
study 


ALEPH 

i 


L3 

m 

prelim. 


OPAL 

|15| 


CDF 
|16| 


D0 



VV+ 

vv- 


4.7 
5.0 


6.4 
7.1 


3.8 
5.0 


4.1 
5.7 


3.5 
5.2 


4.9 
6.1 


AA+ 

AA- 


2.6 
3.7 


7.2 
7.9 


5.6 
3.5 


6.3 
3.8 


3.8 
4.8 


4.7 
5.5 


X1+ 
XI- 


1.8 
2.8 


— 


— 


— 


— 


3.9 
4.5 


X3+ 
X3- 


2.9 
2.8 


6.7 
7.4 


4.0 
3.4 


4.4 
3.8 


— 


4.2 
5.1 


X4+ 
X4- 


4.0 
4.3 


2.9 
4.5 


2.9 
4.8 


3.1 
5.5 




3.9 
4.4 


X6+ 
X6- 


2.8 
1.7 










4.0 
4.3 


U3+ 

U3- 


3.7 

3.5 




6.1 
4.9 


4.1 
5.8 






U4+ 
U4- 


4.4 
4.8 




2.1 
2.9 


2.3 
3.2 







° No comparison is made for scenarios for which only ZEUS sets limits: X2, X5, Ul, U2, U5, U6. 

Table 4: Lower A limits at 95% C.L. from this study compared to equivalent 
results from other experiments. The results of the L3 collaboration are preliminary. 
The XI, X3/U3, X4/U4, and X6 scenarios are denoted LL-LR, LL+RR or VO, 
LR+RL or AO, and RL-RR, respectively, by the LEP and Tevatron experiments. 
Note that CDF also provides limits for fifiqq CI (not shown here) which can be 
combined with the eeqq limits if one assumes lepton flavor universality. 
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CI 


Ao 


Ai 


A2 


A3 


A4 


scenario 




As 


Ae 


A7 


As 


W 


0.219960 


-0.317192+2 


0.102541+4 


0.559212+4 


-0.157739+5 






-0.171440+6 


0.296531+6 


0.381829+7 


0.604081+7 


VA 


0.245590 


0.137604+1 


-0.977350+2 


-0.345909+3 


0.131860+5 






0.684982+4 


-0.160571+6 


-0.532035+5 


0.948424+6 


AA 


0.241619 


0.136626+2 


0.255357+3 


-0.361006+4 


0.477991+4 






0.549759+5 


-0.700729+5 


-0.431092+6 


0.831623+6 


XI 


0.206150 


0.700878+1 


0.510132+2 


-0.102610+4 


0.197029+4 






0.829875+4 


-0.124916+5 


-0.338700+5 


0.571381+5 


X2 


0.198413 


-0.113778+2 


0.181979+3 


0.126775+4 


0.163454+4 






-0.100860+5 


-0.147490+5 


0.448248+5 


0.807606+5 


X3 


0.082637 


-0.608908+1 


0.109116+3 


0.214284+3 


0.131258+3 






-0.396591+3 


-0.110789+3 


0.692834+3 


0.517102+3 


X4 


0.327244 


-0.213799+2 


0.470532+3 


0.356225+4 


0.558752+3 






-0.519826+5 


-0.454742+5 


0.362723+6 


0.623043+6 


X5 


0.258595 


-0.142310+2 


0.218532+3 


0.147937+4 


0.816187+3 






-0.134481+5 


-0.960868+4 


0.720345+5 


0.103106+6 


X6 


0.143683 


-0.486270+1 


0.280939+2 


0.843126+3 


0.235808+4 






-0.623902+4 


-0.141521+5 


0.234297+5 


0.505643+5 


Ul 


0.277925 


0.822634+1 


0.581735+2 


-0.708590+3 


0.828391+3 






0.476597+4 


-0.472480+4 


-0.158164+5 


0.208190+5 


U2 


0.086540 


-0.100733+2 


0.350929+3 


0.100296+4 


0.118597+3 






-0.122795+5 


-0.957897+4 


0.960528+5 


0.159882+6 


U3 


0.016648 


-0.397315+1 


0.241440+3 


0.103296+3 


-0.139481+3 






-0.115302+2 


0.154464+3 


-0.851289+2 


0.437650+3 


U4 


0.145224 


-0.224785+2 


0.824204+3 


0.231869+4 


-0.154694+5 






-0.237649+5 


0.454907+6 


0.537245+5 


-0.486230+7 


U5 


0.108391 


-0.149704+2 


0.498695+3 


0.981530+3 


-0.449846+4 






-0.546963+4 


0.877442+5 


0.215059+3 


-0.652759+6 


U6 


0.174813 


-0.417751 + 1 


0.981713+1 


0.462709+3 


0.134709+4 






-0.284097+4 


-0.650313+4 


0.835725+4 


0.172630+5 



Table 5: Pararneterizations of the functions L(e/A^) — Lmin resulting from the 
unhinned method. For each CI scenario, the nine coefficients Aj define a polynomial 
X^i=o Ai(e/A^)* which has been fitted to L(e/A^) — Lmin in the range where L(e/A^) — 
Lmin < 18. The notation x+n is used as a shorthand for x ■ 10". 
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2.5 3 3.5 4 4.5 2.5 3 3.5 4 4.5 




2.5 3 3.5 4 4.5 2.5 3 3.5 4 4.5 

log,o(Q2[GeV2]) log^oiO'lGeV^]) 

Figure 2: Examples of the relative influence of a CI on the NC DIS cross section 
(P(7/{dx dQ"^) at various fixed values of x. The ratio (SM+CI)/SM of differential 
cross sections for the scenarios VV (top left), AA (top right), XI (bottom left) and 
X3 (bottom right) is plotted for A — 2 TeV . The symbols are used only to label the 
different curves. 



23 




■0.1 0.1 -0.1 0.1 0.2 0.3 




0.2 0.4 -0.6 -0.4 -0.2 0.2 

8/A^ IJeV^) 8/A^ ITeY^) 

Figure 3: Log-likelihood as a function o/e/A^ as obtained by the two methods 
described in Sect. for the scenarios VV (top left), A A (top right), XI (bottom 
left), and X3 (bottom right). For both methods, the values of L — Lmin ore shown. 
Note that the two methods use different statistical approaches and, in addition, dif- 
fer in their treatment of systematic effects. Each figure represents two CI scenarios 
(e = +l, e = -l). 
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ZEUS e*p 1994-97 
VV 

• Data/MC 
statistical errors only 
— 95% C.L limit (e=+1) 
95% C.L limit (e=-1) 
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-0.5 
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3.5 4 4.5 

log,o(Q'[GeV']) 



-0.5 



logio(x) 



Figure 4: Exclusion limits on CI strengths in terms of the corresponding mod- 
ification of the expected distributions of (l^ft) o-nd x (right) in the kinematic 
region used for method 1 (see Sect. and Table^. The dots represent the ratios of 
observed and expected numbers of events and the error bars indicate the statistical 
uncertainties. The histograms show the modification of this ratio for the VV ( top ) 
and A A (bottom) scenarios as obtained with the Aum values given in Table^. The 
solid (dotted) line corresponds to the 95% exclusion limit for e = +1 (e = —1 ). 
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ZEUS e*p 1994-97 
X1 

• Data/MC 
statistical errors only 
— 95% C.L limit (e=+1) 
95% C.L limit (e=-1) 




(0 1-5 




3.5 



4.5 



-0.5 



CM 

o 

■a 



w 

(0 

o 
E 

CM 

o 

■a 




3.5 4 4.5 

log,o(Q'[GeV']) 



-0.5 
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Figure 5: Exclusion limits on CI strengths in terms of the corresponding mod- 
ification of the expected distributions of (l^ft) o-nd x (right) in the kinematic 
region used for method 1 (see Sect. and Table^. The dots represent the ratios of 
observed and expected numbers of events and the error bars indicate the statistical 
uncertainties. The histograms show the modification of this ratio for the XI (top) 
and X3 (bottom) scenarios as obtained with the Aum values given in Table The 
solid (dotted) line corresponds to the 95% exclusion limit for e = +1 (e = —1 ). 
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A;;,in(Tev) 



5.0 
2.6 
3.7 
2.8 
3.1 
2.8 
4.3 
3.3 
1.7 
2.6 
3.9 
3.5 
4.8 
4.2 
1.8 



ZEUS e> 1994-97 
(95% C.L.) 



i9- 



4^ 



A^inCTeV) 



4.7 
2.5 
2.6 
1.8 
3.4 
2.9 
4.0 
3.5 
2.8 
2.0 
4.0 
3.7 
4.4 
4.0 
2.4 



■0.4 -0.3 -0.2 -0.1 



0.1 0.2 0.3 

(TeV^) 



Figure 6: Confidence intervals o/e/A^ at 95% C.L. for all CI scenarios studied 

in this paper (horizontal bars). The filled (open) circles indicate the positions of 
the log-likelihood minima, e/A^, for e = +1 (e = —I)- The numbers at the right 
(left) margin are the lower A limits for e = +1 (e = —I). 
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Figure 7: The junction R defined in eq. ( fi^ j for the VV (top) and the AA 
(bottom) scenarios, for e = +1 (solid lines), and for e = — 1 (dashed lines). The 
graphs are based on the polynomial parameterizations of Table |^. 
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